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ELASTIC  RESPONSE  OF  A  CYLINDER 
CONTAINING  LONGITUDINAL  STIFFENERS 


1.  INTRODUCTION 


Cylindrical  shell  theories  have  been  developed  over  a  number  of  years  and  there  exists  a 
large  quantity  of  these  theories  (reference  1)  to  model  the  response  of  thin  shells  to  free-wave 
propagation  and  various  external  loading  conditions.  Thin  shell  theories  usually  involve  a 
flexural  wave  propagating  in  the  shell’s  axial  direction  and  they  sometimes  include  waves 
propagating  in  the  circumferential  direction.  Fully-elastic  isotropic  (thick)  shell  theory  has  also 
been  developed  (reference  2)  and  this  allows  for  higher-order  wave  motion  in  the  radial 
direction,  which  results  in  significant  displacement  and  stress  differences  in  this  direction. 
Transverscly-isotropic  shell  theory  has  been  formulated  (reference  3)  for  shells  that  ha\e 
different  moduli  in  the  axial  direction  than  the  radial  and  circumferential  direction.  The  reaction 
of  orthotropic  cylinder  response  has  been  derived  (reference  4)  for  shells  that  have  different 
moduli  in  all  three  of  the  shell’s  primary  directions.  The  above  references  all  consider  the  shell 
material  properties  to  be  homogeneous  in  each  specific  direction,  although  for  the  latter 
documents  (references  3  and  4)  these  properties  do  not  necessarily  have  to  be  equal. 

Shells  can  be  reinforced  with  external  or  internal  stiffeners  to  add  strength  and  reduce 
weight.  Circumferential  stiffeners  have  been  added  to  a  thin  shell  energy  model  (reference  5)  to 
determine  the  system's  stop  and  pass  bands  and  this  was  specifically  applied  to  aircraft 
fuselages.  Circumferential  stiffeners  have  been  added  to  a  thin  shell  finite  element  model 
(reference  6)  and  this  resulted  in  identification  of  stop  and  pass  bands  and  the  amount  of 
coupling  between  the  specific  modes.  Free  vibration  of  longitudinally-stiffened  thin  cylindrical 
shells  was  investigated  using  Donnell  shell  equations  (reference  7)  and  the  results  were 
compared  to  an  analytical  solution  generated  using  a  Ritz  solution  method.  Axial  stiffeners  have 
been  added  to  a  thin  shell  energy  method  (reference  8)  to  determine  the  system’s  pass  and  stop 
bands  using  a  variety  of  different  stiffener  shapes.  A  study  of  Timoshenko-Mindlin  type  shell 
with  an  internal,  lengthwise  rib  has  been  undertaken  (reference  9)  to  compute  the  resonances  as  a 
function  of  wave  radius  and  shell  angle.  An  analytical  method  for  longitudinally-stiffened  shells 
to  predict  modal  density  and  radiation  efficiency  for  use  with  statistical  energy  analysis  has  been 
developed  recently  (reference  10).  Note  that  references  5  to  10  arc  all  some  type  of  thin  shell 
theory  and  the  reinforcement  is  external  to  the  cylinder  medium. 

This  report  derives  an  analytical  model  of  a  fully-elastic  cylindrical  shell  that  contains  a 
longitudinal  stiffener  at  some  location  in  the  shell’s  interior.  The  governing  equations  of  the 
shell  are  the  Navier-Cauchy  equations  of  motion  in  cylindrical  coordinates  and  the  governing 
equations  of  the  stiffener  are  the  bar  (wave)  equation  in  the  longitudinal  direction  and  mass 
multiplied  by  acceleration  in  the  radial  and  circumferential  directions.  The  shell  is  divided  into 
two  separate  cylinders  so  that  the  stiffener  lies  on  the  exterior  surface  of  the  inner  cylinder  and 
the  interior  surface  of  the  outer  cylinder.  This  allows  the  displacements  to  be  written  without  the 
stiffener  forces  as  functions  of  the  radial,  circumferential,  and  longitudinal  directions  with 


unknown  wave  propagation  coefficients  on  both  of  the  mediums.  Once  these  are  formulated, 
twelve  boundary  value  equations  are  written  using  these  displacements  or  their  derivatives.  The 
equations  that  couple  the  inner  and  outer  cylinder  stresses  contain  the  forces  of  the  stiffener. 
These  equations  are  solved  using  an  orthogonalization  procedure,  and  this  results  in  the  solution 
of  the  wave  propagation  coefficients.  These  are  inserted  back  into  the  original  displacement 
equations  and  this  produces  known  displacements  and  stresses.  The  model  is  validated  by 
comparing  a  zero  wavenumber  (plane  strain)  excitation  to  a  solution  obtained  using  finite 
element  analysis.  An  example  problem  is  formulated  and  analyzed.  The  addition  of  multiple 
stiffeners  into  the  model  is  also  discussed. 


2.  SYSTEM  MODEL 


The  system  model  is  that  of  a  cylinder  containing  a  longitudinal  stiffener,  as  shown  in 
figure  1 .  This  problem  is  analytically  modeled  by  assuming  that  the  cylinder  is  governed  by 
fully-elastic  equations  of  motion  and  the  stiffener  is  modeled  using  the  bar  wave  equation  for  its 
axial  (extensional)  motion  and  line  mass  equations  for  its  non-axial  motion.  The  model  uses  the 
following  assumptions:  (1)  the  cylinder  has  infinite  spatial  extent  in  the  axial  direction;  (2)  the 
displacements  in  the  cylinder  are  linear  and  three-dimensional;  (3)  the  longitudinal  stiffener  is 
continuous  and  it  exerts  forces  in  the  radial,  circumferential,  and  axial  directions  on  the  cylinder, 
and;  (4)  the  stiffener’s  cross-sectional  area  is  sufficiently  small  that  it  can  be  modeled  as  a  point 
force  in  the  radial  and  circumferential  directions. 


Figure  1.  Schematic  of  a  Modeled  System 
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The  cylinder’s  motion  is  governed  by  the  Navier-Cauchy  equations  of  motion  written  in 
vector  form  as 


pV2u(r,0,z,t)+(A  +  //)VV«u(r,0,z,A)=/> 


d2u  (r,0,z,t) 

a ? 


(i) 


where  p  is  density,  A  and  p  are  Lame  constants,  r  is  the  radial  direction,  0 is  the  circumferential 
direction,  z  is  the  axial  direction,  t  is  time,  and  u  is  the  cylindrical  coordinate  displacement 
vector.  Using  previously-developed  techniques  (references  1 1  and  12),  and  separating  the 
domain  into  two  regions  based  on  the  radial  location  of  the  stiffener  b ,  the  radial  displacement 
can  be  written  as 

Uj(r,0,z,t)=  ^t/i/)(/')exp(i/»^)exp(iAz)exp(-ir«rf)  ,  (2) 

m=  -  to 

where  a  <r  <b  is  denoted  region  (j  =)  1  and  the  radial  dependency  of  equation  (2)  is  written  as 
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and  where  b  <  r  <  c  is  denoted  region  2  and  the  radial  dependency  of  equation  (2)  is  written  as 
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The  circumferential  (tangential)  displacement  can  be  written  as 
vj(r,0,z,t)=  ^ V ’ ( r) exp(\m0) exp(ifcz)  exp( -i cot )  , 

w--- * 

where  a  <  r  <  b  the  radial  dependency  of  equation  (5)  is  written  as 
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and  where  b  <  r  <c  the  radial  dependency  of  equation  (5)  is  written  as 
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The  axial  (longitudinal)  displacement  can  be  written  as 


Wj(r,0,z,t)=  exp(iwfl)exp(iAz)exp(-ift#)  , 

OT=-0O 


(8) 


where  a  <  r  <b  the  radial  dependency  of  equation  (8)  is  written  as 
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and  where  b  <  r  <c  the  radial  dependency  of  equation  (8)  is  written  as 
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In  the  above  equations,  a  is  the  inner  radius  of  the  cylinder,  b  is  the  radial  location  of  the 
stiffener,  c  is  the  outer  radius  of  the  cylinder,  a  is  the  modified  wavenumber  associated  with  the 
dilatational  wave,  /?  is  the  modified  wavenumber  associated  with  the  shear  wave,  k  is  the  axial 
wavenumber,  co  is  the  frequency,  Jm  is  an  mth  order  standard  Bessel  function  of  the  first  kind, 
Ym  is  an  mth  order  standard  Bessel  function  of  the  second  kind,  and  A„ ,  Bm ,  Cm ,  Dm ,  Em, 
Fm,  Gm,  H„ ,  K „ ,  Lm ,  Mm  and  N„  are  unknown  wave  propagation  coefficients  whose 
solutions  are  determined  below.  The  modified  dilatational  wavenumber  is  calculated  using 
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a  = 


(11) 


where 


cd  = 


A  +  2/y 


and  the  modified  shear  wavenumber  is  calculated  using 


(12) 


P  = 


where 


(13) 


(14) 


Note  that  imaginary  arguments  in  the  Bessel  functions  are  permissible.  For  the  remainder  of  this 
report,  the  exponentials,  with  respect  to  time,  are  suppressed  in  all  of  the  equations. 

The  solution  to  the  wave  propagation  coefficients  are  found  using  12  boundary  value 
equations  that  are  written  using  the  stress-strain  and  displacement  relations  at  the  locations  of  the 
shell’s  inner  radius  ( r  =  a),  location  of  the  stiffener  (r  =  b),  and  the  outer  radius  of  the  shell  (r  = 
c).  At  the  inner  radius  of  the  shell  (/*  =  a),  the  cylinder  is  modeled  as  a  free  surface  and  the 
normal  stress  in  the  radial  direction  is 


,  ^  ,  ,,  „  ,  du,  (a,9,z)  A  ,  „  . 

vrr  (a,  0,  z)  =  (A  +  2/j) - — - -  +  —!/, 

or  a 

dz 


A  dv\(a,0,z) 

a  do 


and  the  shear  stresses  associated  with  the  radial  direction  are 
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At  the  radial  location  of  the  stiffener  (/-  =  b),  three  stress  balances  are  written  that  equate  the 
forces  in  the  embedded  stiffener  to  the  stresses  in  each  region  of  the  cylinder.  The  normal  stress 
in  the  radial  direction  is 
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and  the  shear  stresses  associated  with  the  radial  direction  are 
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where  fr(b,0,z) ,  f0(b,6,z) ,  and  f:(b,6,z)  are  the  radial,  circumferential,  and  longitudinal 
stresses,  respectively,  that  the  embedded  stiffener  exerts  on  the  shell,  8  is  the  Dirac  delta  function 
and  6o  is  the  angular  offset  of  the  stiffener  with  respect  to  6  =  0 .  At  the  stiffener’s  radial 
location,  three  equations  are  written  that  equate  the  displacement  continuity  from  the  first  region 
to  the  second  region.  These  equations  are 


u{(b,0,z)  =  u2{b,0,z)  , 

(21) 

v](b,6,z)  =  v2(b,0,z)  , 

(22) 

and 


W|  (b, 6 ,  z)  =  w2  ( b ,  6 ,  z)  .  (23) 

At  the  outer  radius  of  the  shell  (r  =  c),  the  cylinder  is  modeled  as  a  free  surface  with  an  applied 
radial  stress.  The  normal  stress  in  the  radial  direction  is 
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and  the  shear  stresses  associated  with  the  radial  direction  arc 
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where  p(0,z)  is  the  applied  radial  stress  at  the  location  r  =  c. 

For  the  embedded  stiffener,  the  internal  forces  in  the  radial  and  circumferential  directions 
are  modeled  as  line  masses  multiplied  by  the  acceleration  in  these  directions.  This  is  written  as 
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where  M,  is  the  mass-per-unit  length  of  the  stiffener.  Note  that  if  the  stiffener  has  significant 
bending  stiffness,  the  mass  terms  in  equations  (27)  and  (28)  can  be  replaced  by  a  dynamic  beam 
model  such  as  the  Bemoulli-Euler  beam  equation  or  the  Timoshenko  beam  equation.  The 
internal  forces  in  the  axial  direction  are  modeled  using  the  bar  wave  equation,  and  this  equation 
is  written  as 
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^  1  (b) cxp(i/n 0)  exp(iAz)  , 


(29) 


where  As  is  the  cross-sectional  area  of  the  stiffener,  p,  is  the  density  of  the  stiffener,  and  £,  is 
Young’s  modulus  of  the  stiffener.  For  the  numerical  examples  presented  here,  the  external  force 
on  the  outside  of  the  shell  is  modeled  as  a  ring  pressure  with  axial  wavenumber  A.  and  is  given 
by 


p(0,z)  =  P0  exp(iAz)  . 


(30) 
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Using  the  above  four  equations,  and  the  functional  form  of  the  displacement  fields,  equations 
(15)  to  (26)  become 


I 


d Ul  (a)  ,  A  ,r(l)/_%  ,  \mA  I/(1) 


(A  +  2//)“""  +  +  - 

dr  a  a 


V™(a)  +  ikXW"\a) 


exp(i mO)  =  0  , 


(31) 


dr  a  a 


exp(iw$)  =  0  , 


(32) 


i  kj.iU"\a)  +  fj 


djjgV) 

dr 


exp(i///$)  =  0  , 


(33) 


(A  +  2^)dUl(b)-  +  ^-U^(b)  +  ^Ki'>(6)  +  i  kAW"\b) 

dr  b  b 


exp(iw^) 


-I 


u + 2//i  + 4  c/lr1  (/>) + ^  V?’  (h) + (M 

dr  n  n 


exp(iw^) 


ru-A/t 


2nb 


"L_ 

X  1  (/>)  exp(  iw  ) 


8(0 -0O) , 


(34) 


dr  b  b 


exp(i  mff) 


n 


-i 

m=  -o 

&rA/  m=+” 


2^6 


dr  6 

_+ao 

^fil)(6)exp(iw6') 


exp(imtf) 


8(0-0 n). 


(35) 


Z 


iAyr(/!»  +  // 


d^°(6) 


-I 


dr 

ikpU{J)(b)  +  p 


exp(i//7$) 
d  W?\b) 

dr 


exp(i///$) 


(A,psco2  -  A,Esk2) 


2  ub 


^^l,(6)exp(iwt9) 


8(0  —  0ft) , 


(36) 
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(37) 


n- 

£t/i"(6)exp(im0)-  ^(7j„2)(/>)exp(iw0)  =  0 


m=  oc 


m=  x 


w=  +  x>  m=+x> 

2>i>)exp(im0)-  X  exp(im<9)  =  0  ,  (38) 


H  =  +  X  W=+X 

X^m'(6)cxp(im^)-  £Vj,2,(/>)exp(iW0)  =  O  ,  (39) 


(A  +  2/y) - - - 


X 


i  7(2)/  x 


i  m  / 


exp(  i/7/79)  =  /^,  ,  (40) 


dr  c  c 


exp(iwtf)  =  0 


(41) 


and 


tf=  +  T< 

I 


\kfjU'£)(c)  + fj 


dJrj2,(c) 

dr 


exp(iw#)  =  0  . 


(42) 


The  presence  of  the  Dirac  delta  functions  in  equations  (34)  to  (36)  complicates  the  form 
of  these  equations.  Converting  the  delta  function  to  its  equivalent  Fourier  series  and  rewriting 
the  stiffener  force  terms  will  allow  the  equations  to  decouple  using  an  orthogonalization 
procedure.  First,  the  delta  function  is  rewritten  using  its  Fourier  expansion  on  0  e  [0,2;r]  as 


1  ff=+QO  I  //=+0O 

6(#-6>„)  =  —  Yexp[m(<9- #„)]=—  Yexp(-iM<9o)exp(i/?#)  . 
2 n  n=  r  2 n 

Next,  the  identities 


(43) 
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m=+x 

|  tl=+0C 

£ < 

>exp(i  md) 

—  Vexp(-i/7#0)exp(w#) 

_2^  „=-« 

m=-x 

W"\b) 

i  /»=+»  m=+x 

—  V  V 
2  n£» 


^*(6) 

V2'(b) 

WlX)(b) 


exp[i(/7  -  m)^0  ]exp(im6^) 


(44) 


are  used  to  convert  the  double-multiplicative  summations  into  double-embedded  summations. 
The  proof  of  equation  (44)  is  shown  in  appendix  A.  These  equations  are  substituted  into  the 
right-hand  terms  of  equations  (34)  to  (36),  and  equations  (31)  to  (42)  are  multiplied  by  the 
exponential  exp(-i/?tf)  and  the  resulting  expressions  are  integrated  over  [0,2;r].  Because  the 
exponential  functions  exp(-i/?#)  and  exp(im^)  are  orthogonal  on  this  interval  when  p  , 
equations  (3 1 )  to  (42)  decouple  into  an  infinite  number  of  individually  /^-indexed  equations  given 
by 


(A +  2  M) 


d g(fl) 

dr 


+-Ui"(a)+'-^V},l'(a)  +  ikAW")(a)  =  0  , 
a  a 


(45) 


dr;n(a) 

'IT 


a  a 


(46) 


\k{jUy(a)  +  /j 


dW"\a) 

dr 


=  0 


(47) 


dr  b  b 

-(A  +  TUr] ~  X£T V?2)(b)  ~  i kAttf](b) 

dr  b  b 

=  Z  exP[H«  -  p)0 0  V » ’  (b) ,  (48) 

2.7W  „=  _ 
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d  V"\b) 

T~ 


-^V"\b)+^u"\b) 


d  Vf'(h) 
dr 


+Eyf\b)-mulp2\b) 


(49) 

dr  dr 

( Aspsa r  -  AsEsk2)  lu/(1)/i.x 

=  -  -  2^exPt1("  P)0«W»  (b), 

2  nh 

(50) 

U\]\b)-U{2\b)  =  0  , 

(51) 

s 

1 

tJ 

S- 

ii 

o 

(52) 

o 

II 

C* 

1 

■s 

(53) 

(A  +  2//)d^  >(C)  +—Uip){c)  +  ^  V«'(c)  +  i  kAWf'ic)  =  S0ph , 
dr  c  c 

(54) 

^  d^!,(c)  _  fi  pm(c)+  VU  y«.(£.)  =  0  _ 
dr  c  c 

(55) 

and 


i k/jUp  ( c)  +  p 


dr 


=  0 


(56) 


Note  that  6„r  in  equation  (54)  is  the  Kronecker  delta  function. 


The  left-hand  side  of  equations  (45)  to  (56)  are  /^-indexed  and  these  terms  correspond  to 
the  dynamics  of  the  shell.  The  right-hand  side  the  equations  (48)  to  (50)  are  /7-summed  and  these 
terms  correspond  to  the  stiffener  acting  on  the  shell.  The  functional  forms  of  the  displacements 


listed  in  equations  (3),  (4),  (6),  (7),  (9)  and  (10)  are  now  inserted  into  equations  (45)  to  (56),  and 
this  results  in  an  the  algebraic  matrix  equation  for  each  individual  /?-index  written  as 


[A(/>)]{x,}  =  80p{p}  +  X[F(«,/>)]{x„}  ,  (57) 

n=-oo 

where  [A(/?)]  is  the  12  x  12  matrix  that  models  the  dynamic  response  of  the  shell,  {x^,}  is  the 

12  x  1  vector  of  unknown  wave  propagation  coefficients,  [F(/?,  /?)]  is  the  12  x  12  matrix  that 
models  the  interaction  of  the  forces  in  the  shell  and  the  forces  in  the  stiffener,  and  {p}  is  the  12  x 
1  vector  that  models  the  applied  external  load.  The  entries  of  the  matrices  and  vectors  in 
equation  (57)  are  listed  in  appendix  B.  Equation  (57)  is  now  written  for  all  values  of  the  /j-index, 
and  this  results  in 

A  x  =  F  x  +  p  ,  (58) 

where  A  is  a  block  diagonal  matrix  equal  to 


A  =  ••• 


A( — 1)  0  0 

0  A(0)  0  ••• 

0  0  A(l) 


(59) 


with  0  a  12  x  12  matrix  whose  entries  are  all  zero,  F  is  a  rank-deficient,  block-partitioned  matrix 
equal  to 


F  = 


F( — 1, — 1)  F(0,-1)  F(l-l) 

F( — 1,0)  F(0,0)  F(1,0) 

F(— 1,1)  F(0,1)  F(l,l) 


(60) 


p  is  the  system  excitation  vector  that  models  a  normal  ring  load  and  is  equal  to 

P  =  [-  0T  {p)T  0T  .]T,  (61) 

with  0  as  a  12  x  1  vector  whose  entries  are  all  zero,  and  x  is  the  vector  of  unknown  wave 
propagation  coefficients  and  is  equal  to 

*  =  [-  I*i!t  {x„}t  {x,}t  -]t,  (62) 
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where 


x„={4,  Bn  C„  D„  E„  F„  Gn  H„  Kn  Ln  M„  Nn}r  .  (63) 

Note  that  from  equation  (58),  the  addition  of  the  stiffener  force  terms  couples  all  of  the 
displacement  modes  together  and  they  do  not  respond  independently,  even  though  the  equations 
themselves  are  decoupled.  The  solution  to  the  wave  propagation  coefficients  in  equation  (62)  is 
now  found  by  truncating  the  matrices  and  vectors  in  equation  (58)  to  a  finite  number  of  terms 
and  solving 

x  =  [A  -  F]  p  .  (64) 

Once  the  wave  propagation  coefficients  arc  known,  the  displacement  in  all  three  directions  and 
the  stress  distribution  in  the  cylinder  can  be  calculated. 


3.  MODEL  VALIDATION 


The  model  that  was  developed  in  section  2  can  be  validated  for  zero  wavenumber 
response  (plane  strain)  by  comparing  the  analytical  results  to  model  results  calculated  using  a 
finite  element  program.  This  comparison  used  a  shell  where  the  inner  radius  is  0.0825  m,  the 
outer  radius  is  0.0953  m,  Young’s  modulus  is  1.55  x  10  N  m  \  Poisson’  ratio  is  0.45,  density  is 
1250  kg  m  ,  mass-per-unit  length  of  the  stiffener  is  0.75  kg  m  ',  the  angular  offset  of  the 
stiffener  is  n/b  radians,  radial  location  of  the  stiffener  is  0.0889  m,  and  the  frequency  of 
excitation  is  100  Hz.  For  the  zero  wavenumber  case,  the  displacements  in  the  axial  direction  are 
all  zero  and  the  stiffener’s  axial  properties  do  not  enter  into  the  calculations.  The  finite  element 
model  was  run  using  Abaqus  version  6.10  finite  element  program.  The  model  had  6000  plane 
strain  elements,  one  mass  element,  6600  nodes,  and  was  analyzed  using  the  two-dimensional 
standard  steady  state  dynamics  solver. 

Figure  2  shows  a  plot  of  the  transfer  function  of  displacement  divided  by  pressure  versus 
angle  at  the  radial  location  of  0.0889  m.  The  solid  line  is  the  radial  displacement  of  the 
analytical  model  (equation  (2)),  the  dashed  line  is  the  circumferential  displacement  of  the 
analytical  model  (equation  (5)),  the  solid  markers  are  the  radial  displacement  of  the  finite 
element  model,  and  the  circular  markers  are  the  circumferential  displacement  of  the  finite 
element  model.  The  analytical  model  was  calculated  using  27  terms  ( - 1 3  <  ni  <  1 3 ),  which 
produced  a  324  x  324  system  matrix.  For  this  validation  problem,  there  is  broad  based 
agreement  between  the  analytical  model  and  the  finite  element  model  displacements  at  all 
locations  on  the  cylinder. 
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<r,0t)/Po  (m  Pa  )  <r,d/)IPQ  (m  Pa“) 


-8 

x  10 


Angle  9 (radians) 


Figure  2.  Transfer  Function  of  Radial  Displacement  Divided  by  Pressure  vs  Angle  (top)  and 
Circumferential  Displacement  Divided  by  Pressure  vs  Angle  (bottom) 


Note:  In  figure  2,  the  solid  line  is  the  analytical  model,  the  square 
markers  are  the  finite  element  results,  and  the  dashed  line  is  the 
response  of  the  cylinder  without  the  stiffener. 
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4.  THREE-DIMENSIONAL  ELASTIC  EXAMPLE 


A  three-dimensional  example  problem  is  now  formulated  and  discussed.  The  same 
parameters  used  for  the  shell  in  the  model  validation  (section  3)  are  used  to  construct  this 
numerical  example.  The  embedded  stiffener  used  in  this  model  is  Kevlar  grade  2d  that  has  a 
density  of  1440  kg  m  ,  a  Young’s  modulus  of  83  x  10  N  m  ,  and  a  mass-per-unit  length  of 
0.00285  kg  m  .  The  stiffener  has  a  diameter  of  0.00159  m  and  is  in  the  same  physical  location 
as  the  stiffener  in  the  model  in  section  3.  Figure  3  shows  a  plot  of  the  magnitude  of  the  transfer 
functions  of  displacement  divided  by  pressure  versus  wavenumber  at  the  location  r  =  0.0889  m, 
9=  tc/6  radians,  and  z  =  0  m.  The  plot’s  scale  is  in  decibels  and  the  units  are  m  Pa'1.  The  top 
plot  is  the  radial  displacement,  the  middle  plot  is  the  circumferential  displacement,  and  the 
bottom  plot  is  the  longitudinal  displacement.  The  solid  lines  are  the  model  results  with  the 
stiffener  and  the  dashed  lines  are  the  model  results  without  the  stiffener.  The  circumferential 
displacement  for  the  system  without  the  stiffener  is  zero  and  thus  is  not  depicted  on  the  middle 
plot.  The  analytical  model  with  the  stiffeners  was  calculated  using  27  terms  and  the  analytical 
model  without  the  stiffeners  was  calculated  using  a  single  ( n  =  0)  term,  as  the  ring  load  produces 
an  angularly-symmetric  response  for  the  problem  without  the  stiffener  and  higher  order  terms  do 
not  enter  into  the  analysis.  Note  in  the  figures  that  the  resonance  for  the  ring  mode  has  decreased 
in  wavenumber  due  to  the  effects  of  the  stiffener.  Additionally,  the  next  two  modes  of  the 
stiffened  structure  are  visible  as  the  response  is  no  longer  angularly  symmetric. 

Figure  4  is  a  plot  of  the  radial  displacement  divided  by  pressure  of  the  stiffened  system 
versus  wavenumber  and  frequency.  Figure  5  is  a  plot  of  the  circumferential  displacement 
divided  by  pressure  versus  wavenumber  and  frequency.  Figure  6  is  a  plot  of  longitudinal 
pressure  divided  by  pressure  versus  wavenumber  and  frequency.  In  figures  4  to  6,  the  location 
on  the  cylinder  is  r  =  0.0889  m,  9=  n/6  radians,  and  z  =  0  m;  the  scale  of  the  plot  is  in  decibels; 
and  the  units  are  m  Pa  '.  For  this  system,  the  response  is  predominantly  a  combination  of  the  n  = 
0  and  n  =  1  cylinder  modes.  This  was  deduced  by  comparing  the  n  =  0  and  n  =  1  modes  of  the 
unstiffened  system  to  these  figures.  The  wave  speeds  of  the  individual  waves,  however,  have 
increased  with  the  addition  of  the  stiffener.  Note  that  there  is  some  slight  numerical  instabilities 
at  extremely  low  frequencies  and  high  wavenumbers.  This  is  a  location  (in  wavenumber  and 
frequency)  that  is  not  typically  analyzed  because  there  is  no  free-wave  response  in  this  area. 

Three  other  model  attributes  are  briefly  discussed.  First,  in  the  event  that  the  cylinder 
contains  Q  stiffeners  at  the  radial  location  r  =  b  and  angular  location  9q ,  the  right-hand  side  of 
equations  (48)  to  (50)  are  rewritten  as 


&  arM, 

h. i7ib 


^exp[i(«  -  p)Oq]Ull)(b) 


(65) 


&  eo2Mq 

fa  2/ib 


£exp[i(w-/>)0JF(I<,)(&) 

T> 


(66) 
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and 


(Aqpveo2  -  AgEgk^) 

2rt> 


n=+2 o 

^exp[i(«  -  p)dq]W{n'\b)  , 

» 


(67) 


where  the  subscript  q  corresponds  to  the  </th  stiffener.  These  equations  allow  additional 
stiffeners  at  r  =  b  to  be  incorporated  into  the  model.  Second,  if  the  stiffeners  are  at  a  radial 
location  r  *  b ,  then  the  solution  of  the  cylinder  in  the  radial  direction  has  to  be  further 
subdivided  at  the  radial  location  of  each  stiffener  and  six  new  equations  have  to  be  added  to  the 
stress  and  continuity  equations  for  each  stiffener  when  r  *  a  or  r  *  c .  These  are  a  radial, 
circumferential,  and  longitudinal  stress  equation  and  a  radial,  circumferential,  and  longitudinal 
displacement  equation.  If  the  stiffener  is  at  the  cylinder  boundary,  i.e.  r  =  a  or  r  =  c,  then  only 
the  stress  equations  are  needed  to  incorporate  the  effects  into  the  model.  Third,  if  the  external 
load  is  some  function  other  than  a  constant  ring  load,  then  this  function  has  to  be  replaced  by  its 
Fourier  series,  then  expanded  and  orthogonalized  in  a  similar  manner  as  to  the  displacement 
fields.  This  process  is  well  known. 
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Figure  3.  Transfer  Function  of  Displacement  Divided  by  Pressure  v.s  Wavenumber 


Note:  In  figure  3,  the  top  plot  is  the  radial  displacement,  the 
middle  plot  is  the  circumferential  displacement,  and  the  bottom 
plot  is  the  longitudinal  displacement.  The  solid  lines  are  the  model 
results  with  the  stiffener  and  the  dashed  lines  are  the  model  results 
without  the  stiffener. 
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Figure  4.  Radial  Displacement  Divided  by  Pressure  vs  Wavenumber  and  Frequency 


Figure  5.  Circumferential  Displacement  Divided  by  Pressure  vs  Wavenumber  and  Frequency 
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Figure  6.  Longitudinal  Displacement  Divided  by  Pressure  v.v  Wavenumber  and  Frequency 


5.  SUMMARY 


An  clastic  analytical  model  of  a  system  that  consists  of  a  cylinder  containing  an  embedded 
longitudinal  stiffener  has  been  derived.  This  model  was  developed  so  that  the  thick  shells  that 
contain  one  or  more  stiffeners  can  be  accurately  modeled  and  analyzed.  This  new  model  was 
verified  for  the  ease  of  plane  strain  behavior  by  comparison  of  the  results  to  finite  element 
analysis  results.  An  example  problem  was  developed  and  the  behavior  of  the  system  at  various 
frequencies  and  wavenumbers  was  illustrated.  The  inelusion  of  multiple  stiffeners  was 
diseussed. 
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APPENDIX  A 

DOUBLE  SUMMATION  IDENTITIES 


The  proof  of  the  double  summation  identity  given  in  equation  (44)  is  presented  in  this 
appendix.  The  identity  is  critical  in  transforming  the  problem  into  a  form  where  orthogonal 
functions  can  decouple  the  modes.  The  first  term  in  equation  (44)  without  the  2 n  coefficient  is 


m-+x 

^Xm  exp(\mO) 

_w=-x 


X  exp(-i/;0„)exp(i/j<9)  , 


(A- 1) 


where  Xm  is  either  (/^’(/j),  V^(b)  or  W^(b) .  This  can  be  written  as 


77 =+00 

z 


in=+cc 

Y<Xm  exp(i mO) 


exp(-i/;^Q)exp(i/;^)  . 


77= -CO 


777  =  — 00 


(A-2) 


Expanding  the  m  indexed  series  results  in 

/)  =  +  X 

+  X  iexp[i(-l)#]+ A'„exp[i(0)6>]+ A',exp[i(l)#]  +  ---]exp(-i77#o)exp(i/7$)  .  (A-3) 

n-  t: 

Next,  expanding  the  n  indexed  series  yields 

— h[ — t-  X  i  exp[i(-l)#]  +  A'0 exp[i(0)$]  +  X,  cxp[i(l)#]-t — ]cxp[-i(-l)^,]exp[i(-l)<9] 

+  [•••  +  X.,  exp[i(- 1  )0]  +  A',,  exp[i(0)tf]  +  A,  exp[i(l)tf]  +  •  •  •]exp[-i(0)6>„]exp[i(0)6>] 

+  [ — h  X  i  cxp[i(-l)#]  +  A''0exp[i(0)^]  +  X,  exp[i(l)#]  +  •••]exp[-i(l)6>0]cxp[i(l)$]  +  •••  , 


(A-4) 


and  multiplying  through  gives 

•••  +  [•••+  X  ,  exp[-i(-l)<90]exp[i(-2)tf]+  A'0  exp[-i(-l)6(,]exp[i(-l)6>] 

+  A',  exp[-i(- 1  )6>rt  ]  exp[i(0)$]  +  •  •  •] 

+  [•••  +  X  i  exp[-i(0)#(l]exp[i(-l)#]  +  Ar()  exp[-i(0)#„]exp[i(0)t9] 

+  A',  exp[-i(0)(?„]exp[i(l)#]  ■+ — ] 

+  [•••  +  X  i  exp[-i(l)#0  ] exp[i(0)$]  +  X„  cxp[-i(l)tf0  ] exp[i(l)6>] 

+  A',  exp[-i(l)0o]exp[i(2)<9]  +  •••]  +  •••  .  (A -5) 


A-l 


Regrouping  equation  (A-5)  on  specific  values  of  Xm  yields  the  equation 


m-_+ao 

•••+  y'X  ,  exp[-i(w-(-l))#0]exp[i/M#]  +  Y^Q  exp[-i(m-0)#0]exp[imtf] 

m=-oc  m=-o o 

+  exp[-i(w- l)^0]exp[iw^]  +  ...  , 

m=-  oo 


which  can  be  rewritten  as 


w=+x  m~+TQ 

Y  Y  exp[i(n- w)0o]exp(iw^)  . 

//=  -x  m--y. 


A-2 


(A -6) 


(A-7) 


APPENDIX  B 

MATRIX  AND  VECTOR  ENTRIES 


The  entries  of  the  matrices  and  vectors  in  equation  (57)  are  listed  below, 
entries  of  [A(/?)]  are 


2  pa 

Q\.\  — - }  p+\(aa)  + 

a 


lMRzR.a'-U  +  2p).^ 


a 


2/70  ,,  /  \ 

«I.2  = - Y p  1(00)4- 


a 


2/jp{p-\)  2  2 

£ - a  (A +  2p)-Ak 


a 


Y ,,(aa)  , 


-2i  ppp  .  ,  n  ,  2\f.tp(p-\)  .  ... 

7/1.3  = - Jp+\(fla)-i - ; - )  p(J3a)  , 

a  a 


-  2i ppp  2i pp(p  - 1) 

t/m  = - Ypti{j3a)  + - ; - -Y P(pa)  , 

a  a~ 


fli.s 


2iM'(g_+ 1)  j  ^ {j3a)  +  2iMf3k J  ,  (aa)  , 
a 


(f3a)  4-  2\p/3k  Yp(aa)  , 
a 


-2\pap>  .  .  .  2\pp(p- 1)  . 

«2.i  = - J^ooJ-f  , 

77  77" 


-2\pocp ..  .  .  2\pp>(p  - 1)  . 

772.2  = - I (077)4-  ; - Yp{aa)  , 

77  Cl  ~ 


7/2.3=-^  „+.(/&)  + 


a 


~2/jp{p~\) 


a~ 


+  mP2 


J/,(/A?)  , 


The  nonzero 

(B-l ) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 

(B-6) 

(B-7) 

(B-8) 

(B-9) 


B-l 


«2.4  =-^-Yp»(J3a)  + 


a 


-2 


a 


Y ,(/&)  , 


(B-10) 


«2.5  = 


-  pkp{\  +  \)-2\f.ik 


a 


}p^(pa)  +  \pPk]p{Pa)  , 


(B-l  1) 


a2.6  — 


-  f.ikp(\  +  i)  -1\pk 


a 


Yp^(/h)  +  ip/lkYp(/h)  , 


(B-12) 


a3J  =-2\kpai  p^(aa)+  ~1~¥.  ]  p(aa)  , 

a 


(B- 13) 


,,  ,  x  2i  ukp  , 

ax 2  =  -2iAy/ffYp+i(ciw)  + - Yp(aw)  , 

a 


(B- 14) 


aXi=-^ip{Pa)  , 
a 


(B-15) 


a3.4  =  -^-Y p(j3a)  , 
a 


(B- 16) 


«3.5  = 


Jp+1(> 3a)-^Jp(/h)  , 
a 


(B-17) 


«3.6  = 


Yp+l(/Ja)-'-^Y „{f3a)  , 
a 


(B-18) 


2  pa  . 

a4.i  =  — —  JP+i  (ab)  + 
b 


2Mp-D_g; u+2 


ip(ab)  , 


(B- 19) 


«4.2=^Yp+1  («/>)  + 

b 


IwtP. _j)  -  a 2  (A  +  2p)  -  Ak2 


Y p{ab)  , 


(B-20) 


B-2 


^4,3  ~ 


^4.4  — 


^4.5  ” 


^4,6 


^4.7  = 


^4,S  “ 


^49  = 


^410 


^4.1  I 


^4,12 


ChA  = 


dy  2  = 


2'PfiP  ,  ,  2\pp(p-\)  ,  //3t, 

— - Jp*i  (/>/’)  + - - - ir(pb) 

b  b~ 


^Y^o»)+SaSfzaY,o») 

o  o 


-  2\f.ik{p  + 1) 


J  „.,(/?/>)  +  2ip/JkJr  (pb) 


-2\pk{p  +  \) 


Yp+l(fJb)+2xpf1kYp(pb) 


~aAA  i 


2 


”#4,3  * 


:  -C74<4  , 


:  ^4  5  , 


^4,6  5 


l\pap  y  2i/jp(p-l),  .  Lx 

— P^J,,  ,(«/>)+  AF  ; — -J,,(or/>)  , 

b  b~ 


3eSEY^(<*)+a^pJlY,(a6)  , 


(B-21 ) 

(B-22) 

(B-23) 

(B-24) 

(B-25) 

(B-26) 

(B-27) 

(B-28) 

(B-29) 

(B-30) 

(B-31) 

(B-32) 


B-3 


#5,3  = 


#5,4  — 


#5,5  = 


#5.6  = 


#5.7  = 


#5,8  = 


#5,9  — 


#5.10 


#5.11 


#5.12 


#6.1  = 


#6,2 


-_2±l 

b 


V.  w+ 


z2±L 

b 


Y  PAPb)  + 


-2fjp(p-\) 

b2 


~  2fjp(p-\) 

b2 


+  M/32 


+  M/J2 


j Pm , 

y Pm , 


-  jukp(  1  +  i)-  2i/z& 


-/z^(l  +  i)-2i/v^ 
b 


Jp+i(/3b)  +  [p/3kJp{/3b)  , 


Yp^(/3b)+\M/3kYp(/Jb) 


— #5.1  9 


#5,2  ? 

— #5.3  > 


— #5.4  9 


:  #5.5  ^ 


=  #5.6  ^ 


-2iA//eyJp,i(or6)+  }  p(ab)  , 

b 


-2lkpaYp.t(ab)  +  ^^-Yp(ab)  , 
b 


(B-33) 

(B-34) 

(B-35) 

(B-36) 

(0-37) 

(B-38) 

(B-39) 

(B-40) 

(B-41) 

(B-42) 

(B-43) 

(B-44) 


B-4 


b 


a,A=-^-\„{pb)  , 
b 


#6.5  - 


//(i-l)pO,+2)  ' 

b  ~ 


b 


#6,6  = 


Mi-Dp(f»  +  2)  i-tl) 

b“ 


Y„W-?Y „G«o  , 


(B-45) 


(B-46) 


(B-47) 


(B-48) 


#6.7  =  ”#6,l  i 

(B-49) 

#6.X  =  #6.2  i 

(B-50) 

#6.9  =  #6.3  i 

(B-5 1 ) 

#6.10  =  #6,4  i 

(B-52) 

CS 

II 

1 

£ 

(B-53) 

#6.12  “  ”#6.6  i 

(B-54) 

5 

cl|  -c 
+ 

+ 

5 

1 

II 

r-’ 

Sj 

(B-55) 

«7.2  =-0'Ypt|(o'/j)  +  —  Y,,(ab)  ,  (B-56) 


B-5 


a^=^ip{j3b)  , 
b 


tf7.4=TYpW  , 

b 


^7.5  =\k]p+\(pb)  , 


a7.6  =  \kYr+](fib)  , 


#7,7  -  #7.!  ? 


#7,8  —  “#7,2  i 


#7,9  —  ”#7,3  ^ 


#7,10  ~  ”#7,4  i 


#7.11  —  ”#7,5  - 


#7,12  —  ”#7.6  , 


«8.i  =^-ip(ab)  , 
b 

a*.2  =^Y„(afr)  , 
b 


a,,3=pjp+i(Pb)-^ip(pb)  , 

b 


(B-57) 

(B-58) 

(B-59) 

(B-60) 

(B-61) 

(B-62) 

(B-63) 

(B-64) 

(B-65) 

(B-66) 

(B-67) 

(B-68) 

(B-69) 


B-6 


(B-70) 


a,A=p\pA(Pb)-^p(Pb)  , 
h 

a^=\kip+[{Pb)  ,  (B-71) 

«w=i*Y,#)  ,  (B-72) 

«K.7  =  -«s.i  ,  (B-73) 

«K.S  =  -«s.2  ,  (B-74) 

aK,9  =  — c/js.3  ,  (B-75) 

«k.i<>  =  -chj  ,  (B-76) 

tfx.n  =  -«s.j  ,  (B-77) 

^8,12  =  — #S.6  1  (B-78) 

i79  |  =  iA  J  ,,(ab)  ,  (B-79) 

a9  2  =  iA  Yp  (ab)  ,  (B-80) 

«9.5  APb)-/np(/]b)  ,  (B-81) 

b 

a,, 6  =——p-^-YP>\(fib)-pYp(fib)  ,  (B-82) 

n 


B-7 


#9.7  “  #9,1  , 


(B-83) 


#9,8  — 


#9,!  I 


#9,12 


#10.7 


#10.8 


#10,9 


#10.10 


#10,1  I 


#10.12 


#11,7 


#11,8 


— #9.2 


:  #9,5  > 


:  #9,6  > 


2  jja 


ip,i(ac)  + 


?MP.7)\-a'-{A  +  2/l)-u* 


iJac)  , 


2  pa 


Y  p+,(ac)  + 


2ME_V_a i(A  +  2v)-M2  Y „(ac)  , 


c  c 


2i^Y„,(/fc)+^£<£^y,(A) , 

c  c* 


-  2\pk(p  + 1) 


J/,.,(/?c)  +  2i/y/3tJ/,(ySc)  , 


-  2\pk(p  + 1) 


YpMPc)  +  2ippkYp(/lc)  , 


2iwj,.,(0c)+2w(fzi)jp(ac) , 


^Y,.,<ac)+*Hfc2>Yl,<«c)  , 

c  c‘ 


(B-84) 

(B-85) 

(B-86) 

(B-87) 

(B-88) 

(B-89) 

(B-90) 

CB-91) 

(B-92) 

(B-93) 

(B-94) 


B-8 


^1 1.9  “ 


(B-95) 


J  „♦.  (>&■)  + 


-2pp(p-\) 


*  rifle)  * 


an,o=~^-Yp+](/k)  + 


-2w(p-\) 


+  pP~ 


Yr(fk)  , 


a\\.u  = 


-  pkp(\  +  \)-2\pk 


ip+x(pc)  +  \ppkip(Pc)  , 


1.12  — 


pkp(  1  +  i)-2i/z& 


Yp+1  (/Xc)  +  \ppk  Yp((k)  , 


t7i2.7  =  -2\kpa  J^i  (ac)  + 


2i  fjkp 
c 


J P(ac)  , 


^I2.K 


-2\kpa  Y;,,|  (ere)  +  —  Y/i(q'c)  , 

c 


atu.-^-Y'Vc)  , 

c 


^1 2,1 1  — 


and 


^12.12  — 


//(i-l)p(jy  +  2) 


+  p{p2-kl) 


\p^fk)-^\r(Pc) 


The  nonzero  entries  of  [F(/j,  p)]  are 


(B-96) 


(B-97) 


(B-98) 


(B-99) 


(B-100) 


(B- 101) 


(B-102) 


(B- 103) 


(B-104) 


B-9 


co2M 


Ai  =  T~r Lexp[i(w-p)^0] 
2  rcb 


-a  J„+l («/>)  +  -  J  „(«/>) 
b 


(B-105) 


(o2M 


A:  =  s  exp[i(/?-p)6>0] 
Am 


-aY„+l(ab)  +  -Yn(ab) 
b 


(B-106) 


co  M,  ....  .  _  , 

Aj  =-z-7-exp[i(/i-p)^,] 

2m 


b 


(B-107) 


A.4  =  —  exp[i(/7  -  /?)6>0] 
2m 


^Y„(/fc) 

b 


(B-108) 


As  =“T~r"exp[i(/7 - p)0o]  [i*  hAPbj\ . 


2^) 


(B-109) 


A.6  =^^exp[i(«-p)0o][i*  Yn+1(/%)]  , 


27ih 


(B-110) 


Ai  =£^exp[i(«-p)^o] 

A7W 


XA-)„(ab) 

b 


(B-lll) 


r  co~Ms  r.,  x  . 

A  2  =  — —  exp[i(K-/7)fl0] 

Am 


~  Y„(ab) 
b 


(B-112) 


y.  CO  M,  ... 

A  J  =  —  exp[i(/7-p)6>o] 
2m 


Pho(pb)—-J„(Pb) 

b 


(B-113) 


A.  = 

2m 


P  Yn+1  (/?/>)-—  Y„(f3b) 
b 


(B-114) 


B-10 


As  =i^exp[i(/7-p)6>0][U  , 

l7W 

2  \  i 

A 6  =^!— -exp[i(«-/>)<9„][i*  Y„tl  (/?/>)]  . 

2  7m 


A,  i 


/l*  pv  acT  —  /4j  A  ‘ 
2/rf?~ 


exp[i(«-/?)(9()][i^  J  „(«/>)]  , 


A,2 


Aspsa> 2  -  AsEsk 2 
2nb 


exp[i(/;-/?)<90][i£  Y„  (<ar^>)]  , 


.A.5  - 


A,ppo2  -  AsE,k 2 
2nb 


exp[i(«-/?)0o] 


(1  -i)» 


Jo+,09A)-/?Jn(/%) 


and 


As  — 


A^psor  -  A^EJr 

2nb 


exp[i(«-p)<90] 


(l-i)» 


Y„,\(/Jh)-f}Yn(/]b) 


The  nonzero  entry  of  p  is 
Pm  =  Po  ■ 


(B-115) 


(B-l  16) 


(B-l  17) 


(B-l  18) 


(B-l 19) 


(B-120) 


(B- 121) 


B-l  1  (B-l 2  blank) 


APPENDIX  C 
LIST  OF  SYMBOLS 


A 

Mp) 

Am 

A, t 

a 

a 

Bm 

b 

P 

Cm 

c 

Cj 

cs 

Dm 

Em 

Es 

F 

F  (n,p) 

Fm 

fr(b,0,z) 

.fo(b,  0,z) 

f:(b,0,z) 

Gm 

Hm 

Km 

k 

Lm 

A 

Mm 

Ms 

m 

P 

Nm 

P 

P 

Po 

P 


Dynamic  shell  matrix 
Dynamic  shell  sub-matrix 

Wave  propagation  coefficient 
Cross  sectional  area  of  stiffener 
Inner  radius  of  cylinder 

Modified  wavenumber  associated  with  the  dilatational  wave  of  the  cylinder 
Wave  propagation  coefficient 
Radial  location  of  stiffener 

Modified  wavenumber  associated  with  the  shear  wave  of  the  cylinder 

Wave  propagation  coefficient 

Outer  radius  of  cylinder 

Dilatational  wave  speed  of  cylinder 

Shear  wave  speed  of  cylinder 

Wave  propagation  coefficient 

Wave  propagation  coefficient 

Young’s  modulus  of  the  stiffener 

Stiffener  force  matrix 
Stiffener  force  sub-matrix 

Wave  propagation  coefficient 
Radial  stress  of  the  stiffener 
Circumferential  stress  of  the  stiffener 
Longitudinal  stress  of  the  stiffener 
Wave  propagation  coefficient 
Wave  propagation  coefficient 
Wave  propagation  coefficient 
Longitudinal  wavenumber 
Wave  propagation  coefficient 
First  Lame  constant  of  cylinder 
Wave  propagation  coefficient 
Mass  per  unit  length  of  the  stiffener 
Cylindrical  mode  number 
Second  Lame  constant  of  cylinder 

Wave  propagation  coefficient 
Load  sub-vector 
Load  vector 

Magnitude  of  applied  pressure 
Equation  index 


C-l 


p(0,z) 

r 

P 


Ps 

t 

6 

do 

T,r(a,0,z) 

T,g(a,0,z) 

Tr.Xa,0,z) 

uXr,0,z,t) 

vj(r,0,z,l) 

Wj(r,0,z,t) 

co 

x 


Applied  pressure  at  the  exterior  of  the  cylinder 
Radial  location  on  the  cylinder 
Cylinder  density 

Stiffener  density 
Time 

Angular  location  on  the  cylinder 
Angular  position  of  stiffener 
Normal  stress  in  radial  direction 
Shear  stress  in  radial-circumferential  direction 
Shear  stress  in  radial-longitudinal  direction 
Radial  displacement  of  the yth  layer 
Circumferential  displacement  of  the  jth  layer 
Longitudinal  displacement  of  the  jth  layer 

Frequency 
Coefficient  vector 
Coefficient  sub-vector 
Longitudinal  location  on  the  cylinder 


C-2 
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